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Abstract 

This paper investigates three open problems in random beamforming based communication systems: 
the scheduling policy with heterogeneous users, the closed form sum rate, and the randomness of 
multiuser diversity with selective feedback. By employing the cumulative distribution function based 
scheduling policy, we guarantee fairness among users as well as obtain multiuser diversity gain in 
the heterogeneous scenario. Under this scheduling framework, the individual sum rate, namely the 
average rate for a given user multiplied by the number of users, is of interest and analyzed under 
different feedback schemes. Firstly, under the full feedback scheme, we derive the closed form individual 
sum rate by employing a decomposition of the probability density function of the selected user's 
signal-to-interference-plus-noise ratio. This technique is employed to further obtain a closed form rate 
approximation with selective feedback in the spatial dimension. The analysis is also extended to random 
beamforming in a wideband OFDMA system with additional selective feedback in the spectral dimension 
wherein only the best beams for the best-L resource blocks are fed back. We utilize extreme value theory 
to examine the randomness of multiuser diversity incurred by selective feedback. Finally, by leveraging 
the tail equivalence method, the multiplicative effect of selective feedback and random observations is 
observed to establish the individual rate scaling. 
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I. Introduction 

In multi- antenna downlink systems, transmission strategies which require less feedback re- 
sources [|T)-[H) to fully utilize multiuser diversity [|5], [J6J , but with asymptotic sum capacity 
comparable to dirty paper coding [T71— iTTTI . are favored. The idea of random beamforming [|T2l . 
which satisfies the two aforementioned features has drawn much interest in recent years [fT3l - 
lfl9l . In the basic random beamforming strategy suggested in [12], the transmitter with M transmit 
antennas generates M random orthonormal beams and requires each user to feed back the SINR 
experienced by them for each beam. Then the transmitter schedules users for transmission that 
currently have the best channel for each random beam. Despite the considerable literature on 
this topic, there are three existing open problems: 

1) How to address heterogeneous users with diverse large scale channel effects and the impact 
on scheduling policy? 

n 

2) What is the closed form sum rate by exaclu performance analysis? 

3) What is the effect of selective feedback, both spatial and spectral, on the randomness of 
multiuser diversity? 

The first problem is related to a practical downlink system setting with asymmetrically located 
users having heterogeneous large scale channel effects. This near-far effect was first treated in 
lfl2l by observing that the system becomes interference dominated when M is large enough. In 
the large M setting, the authors prove that users are asymptotically equiprobable to be scheduled. 
However, when M is finite and not increasing simultaneously with the number of users, the 
greedy scheduling policy employed in Ifl2l can not maintain fairness among users. Also, if a 
round robin scheduling policy was utilized, fairness can be guaranteed, but no multiuser diversity 
gain could be achieved for capacity growth. Therefore, an alternate scheduling policy is needed 
to maintain fairness while exploiting multiuser diversity at the same time. In this paper, the 
cumulative distribution function (CDF)-based scheduling policy 11201 is leveraged and analyzed 
in the random beamforming framework, wherein the user whose rate for a given beam is high 

'We use the term exact to denote results valid for arbitrary but finite number of users as opposed to asymptotic results. 
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enough but least probable to become higher is selected. Under this scheduling policy, each user 
can be equivalently viewed as competing with other users with the same CDF, thus making the 
study of individual user rate more relevant and interesting than that of the sum rate. In this 
paper, we develop the notion of individual sum rate, which is the individual user rate multiplied 
by the number of users, in order to demonstrate the multiuser diversity gain with user growth 
for a given user. 

The second problem addresses exact system analysis, namely deriving closed form expression 
for the sum rate for arbitrary but finite number of users. Note that even with full feedback, 
wherein each user conveys back the signal-to-interference -plus-noise ratio (SINR) for M beams, 
the closed form sum rate has not been derived. This is partially due to the complicated form of 
SINR and its interplay with multiuser diversity. In this work, the problem is tackled and solved by 
a novel probability density function (PDF) decomposition [|2T| which decomposes and interprets 
the selected user's SINR. In fl2TJl, the homogeneous setting is considered and in this paper, the 
technique is extended to the heterogeneous user setting and the closed form individual sum rate 
is derived. The closed form result under full feedback helps in evaluating the system performance 
and acts as the building block for exact analysis with selective feedback. 

The third problem is concerned with standard selective feedback in the spatial dimension, 
wherein each user feeds back the SINR for the best beam among the M beams and the corre- 
sponding beam index. This selective feedback is fundamentally different than full feedback in 
two aspects. The first difference is the two-stage maximization with the first stage maximization 
carried out by each user for feedback selection and the second stage maximization carried out 
by the scheduler to perform user selection. Since the best beam is selected by each user, the first 
stage maximization is over M correlated SINR. This correlation issue has been addressed in [22], 
11231 . and the CDF for the selected SINR at the user side is derived. In this paper, we propose 
an approximation for the CDF and utilize it to derive closed form rate approximation. The other 
fundamental difference is the number of the SINR values that the scheduler has to maximize 
over for each beam. This number is fixed and equals the number of users in the full feedback 
case. However, with selective feedback, it becomes a random quantity. In other words, selective 
feedback results in a random effect on the multiuser diversity. This effect was first observed 
in ll22l . In this paper, we investigate the randomness of multiuser diversity by extremes over 
random samples and provide a rigorous argument on the rate scaling. 
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The third problem is further extended to include spectral selectivity by examining a wideband 
OFDMA system, which groups the subcarriers into resource blocks ll24ll to form the basic 
scheduling and feedback unit. In order to save feedback resource while not significantly degrading 
the system performance, additional selective feedback in the spectral dimension is necessary. 
The effect of random beamforming in a wideband system is examined in ll25l by extensive 
simulations, and further studied from a utility function perspective with the proportional-fair 
scheduler in ll26l . In ll27l . analytical results on the asymptotic cluster size is provided. Apart 
from the thresholding-based partial feedback strategy [|28l . the best-L selective feedback strategy 
E9l is appealing and utilized in practical systems such as LTE 11301 . In this paper, we employ 
the best-L selective feedback strategy to investigate random beamforming and the effect of 
spectral dimension selective feedback, which calls for an additional maximization stage at the 
user side to perform feedback selection. In this feedback strategy, only the best beams from 
the best L resource blocks along with the beam and resource block index are fed back from 
each user. In this paper, we first derive a closed form rate approximation with exact analysis, 
i.e., valid for arbitrary but finite number of users. Then, the influence of the additional spectral 
dimension selective feedback on the type of convergence is investigated with the technique of tail 
equivalence. Moreover, the multiplicative effect of selective feedback and random observations 
is observed to establish the rate scaling. 

To summarize, the main contributions of this paper are threefold: the utilization of CDF-based 
scheduling policy to address heterogeneous users with random beamforming, the obtained closed 
form rate results with different selective feedback assumptions, and the asymptotic analysis on 
the randomness of multiuser diversity incurred by selective feedback. These three contributions 
analytically examine the raised open problems, and foster further understanding on random 
beamforming with heterogeneous users and selective feedback. The organization of this paper is 
as follows. Section HI] reviews the basic narrowband system model for random beamforming. The 
analysis for the full feedback case is carried out in Section [Till and for the spatial dimension 
selective feedback in Section [IV] Section |V] provides the model for the wideband OFDMA 
with random beamforming, and examines the effect of additional spectral dimension selective 
feedback on rate performance. Finally, Section |VT] concludes the paper. 
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II. System Model 

We consider a multi-antenna narrowband Gaussian downlink channel with K single antenna 
receivers and a transmitter equipped with M antennas. A block fading channel model with coher- 
ence interval T is assumed. The random beamforming strategy employs M random orthonormal 
vectors (f) m G C Mxl for m — 1, . . . , M, where the 0/s are drawn from an isotropic distribution 
independently every T channel uses |[T2||. Denoting s m (t) as the mth transmission symbol at 
time t, the transmitted vector of symbols at time t, represented by s(t) G C Mxl , is given as 

M 

<t) = J2^(t>m(t), t=l,...,T. (1) 

m=l 

Let Uk(t) be the received signal at the kth user, then 

M 

Vk(t) = y/P^tiWmMsmit) + V k (t), (2) 

m=l 

where h k G C Mxl is the complex channel vector which is assumed to be known at the receiver, 
Vf. is the additive white noise, and the elements of and Vk are i.i.d. complex Gaussian with 
zero mean and unit variance CJ\f(0, 1). Note that this channel assumption corresponds to the 
Rayleigh fading assumption for the small scale channel effect. From now on, the time variable 
t will be dropped for notational convenience. The total transmit power is chosen to be 1, i.e., 
E[s^s] = 1, and thus the received signal-to-noise ratio (SNR) of user k is pu- In a practical 
downlink setting, due to different locations of users, the large scale channel effects pk which 
may consist of path loss and shadowing vary across users. From ©, the SINR of the A;th user 
for the mth transmit beam can be computed as 

SINR fcim = — fefe^— , m = l,...,M. (3) 

M/pk + ^m\^\ 2 

Denote Zk, m — SINRfc iTO for notational simplicity. Then for a given beam m, the ^fc, m 's are 
independent across users k but non-identically distributed due to different pk. For a given user 
k, the Zfc, m 's are identically distributed and correlated. Thus the beam index m can be dropped 
in the expression for the PDF, which is computed in [[12] as 

— — X 

^w = (iSw& +I » + M - 1 )" w ' (4) 
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where u(-) is the Heaviside step function. The CDF of Z k is shown in [fT2ll to be 



Fz h (x)= 1^— ) M 0c). (5) 



III. Full Feedback Analysis 

This section is devoted to the analysis for the full feedback case wherein each user feeds 
back the SINR for M beams. Since under full feedback, all the beams are fed back, the order 
statistics for each beam is over K independent random variables. Thus this case is well suited 
for illustration of the scheduling policy and the derivation of the individual sum rate. 

A. Scheduling Policy and Individual Sum Rate 

After receiving the SINR fc m from user k for beam m, the scheduler is ready to perform 
user selection. In a homogeneous setting, selecting the user with the largest SINR for a given 
beam maintains fairness and obtains multiuser diversity gain. This system was analyzed in our 
recent work [|2T). The work is now expanded to the more complex heterogeneous case. In a 
heterogeneous setting, the greedy scheduling policy would be highly unfair for finite M. The 
round robin scheduling policy can maintain scheduling fairness, but no multiuser diversity gain 
can be obtained. The proportional-fair scheduling policy flU, ll26l achieves the system fairness 
in terms of system utility. However, under the scenario of inter-beam interference, the users' 
rates are coupled under the proportional-fair scheduling policy. This coupled effect makes it 
very difficult, if not impossible, to develop further analytical results^. Therefore, to tackle this 
problem it is useful to consider alternate scheduling policies that decouple each user's rate. In 
this paper, we employ the CDF-based scheduling policy 11201 for further analysis. According to 
this policy, the scheduler will utilize the distribution of the received SINR, i.e., Fz k . It is assumed 
that the scheduler perfectly knows the CDF|j, and it performs the following transformation [20]: 

Zk,m = Fz k {Zk,m)- (6) 

2 Note that extensive simulation results have been provided regarding the use of proportional-fair scheduling policy under 
random beamforming in existing literature such as [26 1 . However, the coupled effect of user's rate prevents further analysis and 
it remains an open problem to theoretically understand the system performance of proportional-fair scheduling policy under the 
heterogeneous user setting with inter-beam interference. 

3 This is the only system requirement to conduct the CDF-based scheduling, and the CDF can be obtained by infrequent 
feedback from users and learned by the system. 
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The transformed random variable Z fe m is uniformly distributed ranging from to 1, and in- 
dependent and identically distributed (i.i.d.) across users for a given beam. Denote k* m as the 
random variable representing the selected user for beam m, then 



k* m = max Z k;m , (7) 



where U m denotes the set of users conveying feedback for beam m. In the full feedback case, 
\U m \ = K. After user k* m is selected per ©, the scheduler utilizes the corresponding >m for 
rate matching of the selected user. Let X m be the SINR of the selected user for beam m and 
now consider the sum rate of the system defined as follows, 



R = E 



M 



(8) 



io g 2 (i + x m 

_m=l 

From the aforementioned formulation, the sum rate can be computed in the following procedure 

R - MEk - l lo ^( 1+F zLJ x) ) dxK 

= #E r^ g2 (l + t)d(F z M K = ^f:MK), (9) 
fc=i ^° k=i 

where (a) follows from the sufficient small probability that multiple beams are assigned to the 
same user; (b) follows from the change of variable x = Fz k , m {t), the fair property of the 
CDF-based scheduling, and the following definition for J k {e) with exponent e E N + : 

POO 

J k (e)± / log 2 (l + x)d(F Zk (x)Y. (10) 



With the help of the CDF-based scheduling, each user feels as if the other users have the same 
CDF for scheduling competition [|20l . Therefore, each user's rate is independent of other users 
making it possible to consider or predict individual user's rate by only examining its own CDF. 
It is clear that the scheduling policy is not only fair, but also acknowledges multiuser diversity 
at the same time. If we denote the sum rate as the "macro" level understanding of the system 
performance, then the individual user rate can be seen as the "micro" level understanding of the 
system performance since this performance metric examines the rate for any specific user and 
the sum rate can be directly computed from the individual user rate from all the users. Thus, 
under the CDF-based scheduling policy, each user's rate can be examined separately and this 
property serves as one building block for further analysis with selective feedback. 
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In order to demonstrate the multiuser diversity gain for each individual user, we define the 
individual sum rate Rk for user k which is the individual user rate Rk multiplied by the number 
of users, as follows 

Rk = KR k = MJ k (K). (11) 

The definition of the individual sum rate under the CDF-based scheduling policy makes it 
natural to examine the rate scaling for each user separately, and also provide a "micro" level 
understanding of the sum rate scaling. Compared with the sum rate and the individual user rate 
which can be treated as performance metrics, the notion of individual sum rate can be regarded 
as the analytic metric for further scaling analysis. 

Note that in the homogeneous setting, Jk{e) reduces to J(t) = / °°log 2 (l + x)d(F z (x)) e . It 
is mentioned in previous works that the exact closed form for J(e) is hard to obtain due to the 
coupled effect of SINR and multiuser scheduling. In the sequel, the closed form expression for 
J7fc(e) is obtained which is the key to computing the sum rate given by ©. The main technique 
is employing the following proposed PDF decomposition which readily follows from [|2T|. 

Lemma 1. (PDF Decomposition) d(F Zk (x)Y can be decomposed as 

e-l N • / M(i+i)x \ 

d( F ZtM) - = e g ( e - - (1 J . az> 

With the help of this PDF decomposition, J7fe(e) can be computed in closed form using standard 
integration techniques whose expression is presented in the following theorem. 

Theorem 1. (Closed Form of Jk) 

™-£g( e 7^(^^- 1M4+ H- (13) 

where = J °° dx whose closed form expression is presented in Appendix lAl 

Proof: The proof is given in Appendix [A] ■ 
Remark: A few remarks are in order. Firstly, the analytically useful PDF decomposition 
decouples the effect of multiuser diversity and random beamforming, which facilitates the 
integration. The decomposition is general in that it can be applied to other channel models, 
though in this paper the simple Rayleigh channel model is assumed to obtain the SINR statistics 
in ©. Secondly, the derived closed from results for the individual sum rate and the sum rate 
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only involve finite sums and factorials, which can readily be computed. Moreover, the derived 
Jkif) will be employed as a building block for rate computation in Section [IV] and Section |V] 
with selective feedback. 

B. Individual Scaling Laws 

With homogeneous setting, the asymptotic sum rate scaling is of interest and has been 
established as M\og 2 \og 2 K Ifl2ll given the SINR statistics in ©. It can be easily seen that 
the multiuser diversity gain is linear with respect to the number of transmit antennas. With 
heterogeneous setting employing the CDF-based scheduling, the same technique can be applied 
to obtain the asymptotic scaling for the individual sum rate R k of user k. We now develop the 
notion of individual rate scaling and state the individual scaling laws under full feedback through 
the following theorem. 

Theorem 2. (Individual Scaling Laws Under Full Feedback) 

lim — R * = 1. (14) 

K^oo M log 2 log 2 K 

Remark: It is seen from Theorem |2] that users asymptotically follow the same scaling laws 
in the CDF-based scheduling policy. The large scale channel effect p k is not written explicitly 
in (fT4l) since it is a constant inside the log term. It should briefly be noted that the rate scaling 
only measures the asymptotic trend when K — V oo and thus can not accurately match the exact 
performance for finite regions of K. 

IV. Selective Feedback in the Spatial Dimension 

This section examines selective feedback in the spatial dimension wherein each user only 
conveys the best beam. This standard user side selection requires the handling of correlated 
random variables and the random effect on observations, which are pursued in Section ITV-AI and 
Section ITV-Bl 

A. Individual Sum Rate 

With selective feedback, each user selects and feeds back the largest SINR among M beams. 
As discussed in Section HH the Z^'s are correlated random variables given k. Thus simple 
order statistics result can not be used to characterize the selected SINR at user side. Denote 
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Yk,m*(k) — niaxZj.^ representing the selected SINR for user k with m*(k) as the selected beam 



rn 



index. Then according to the derivation in 11221 . [23], the CDF of Yk >m *(k) is shown to be 



i=i 

where d % (x) = ^^z^i > A(#) = (40*0 nf^(^( x ) — di(x)), and [•] + is the positive part of 
the argument. Note that the distribution does not depend on the selected beam index m*(k) 
due to the identically distributed property across beams and is dropped to simplify notation, 
i.e., FY km t W {x) = Fy h (x). Using a similar procedure to that described in Section HII-A[ after 
receiving feedback, the scheduler performs the transformation for user selection: 

Yk,m*(k) = Fy h (Yk,m*(k)) ■ (16) 

Compared with ©, it is clear that F Yk = F Zk for the full feedback case. Denote k* n as the 
random variable representing the selected user for beam m, then 

k* m = max Kfc >m *(fc), (17) 

Vim 

where U m = {k : m*(k) = m} denotes the set of users conveying feedback for beam m. lA m is 
a set of random size and the probability mass function (PMF) can be shown to be given by 

K\ ( 1 \ T1 / 1 x K ~ T1 



P(W = n)=y [i-jt) , 0<n<K. (18) 

Following the derivation in Section IIII-AL let X m be the selected SINR for beam m at the 
scheduler side, then conditioned on k* m and \U m \ = T\, the conditional CDF of X m can be 
written as F Xm \ kUUm \ =T1 (x) = (F Yk ^ m (x)) T1 . By averaging over the randomness of \U m \, the 
conditional CDF is expressed as 

F x ^Jx) = t (s)" i 1 - JiT" i^-Wl" <19) 
From © and (fTTTl . the individual sum rate of user k is derived asc 

K / ts-\ / 1 \ T i / 1 \ K—Tl [-CO 

^ = M E( Ti )(m) ^-Ji) I log 2 (l + x)^(x)r.. (20) 



4 In this paper, it is assumed that if no user feeds back SINR for a certain beam, that beam would be in scheduling outage 
and would not contribute to rate calculation. 
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Exact CDF 

Negative Association Upper Bound 
Frechet Upper Bound 



5 10 15 

X 

Fig. 1. Comparison of the exact CDF Fy k with the Frechet upper bound and the negative association upper bound for spatial 
dimension selective feedback (M — 4, pu = 10 dB). 

Due to the complicated form of Fy k , the exact closed form expression for (|20T ) is hard to 
obtain. We now aim to provide an approximate expression for the closed form by examining 
the property of F Yk and utilizing the established result in Section IIII-Al Recall that Y k is the 
maximization over M correlated random variables Z k<m , thus alternative approximation for F Yk 
would lead to rate approximation. One simple approach is to use the Frechet upper bound [1311 
for the Z k>m 's. Since the Z k ^ m 's are identically distributed across m, the Frechet upper bound 
yields F Zk . This upper bound is very loose empirically for F Yk . One suitable approach is inspired 
by the conjectured negative associated upper bound proposed in ll22l to deal with the minimum 
mean square error (MMSE) receiver. Our empirical evidence shows that even with single antenna 
receiver, the Z k ^% are negative associated [32], thus the upper bound produced by the negative 
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association property can be utilized to approximate F Yk , namely 

F Yk (x) ~ (F Zk (x)) M . (21) 

Fig. \T\ illustrates the bounds and the empirical CDF Fy k for M = 4, pk = 10 dB. It can be seen 
that the proposed upper bound in (1211 approximates the exact one in (fl"5l ) well, especially when 
the SINR is large. By using the CDF approximation, the individual sum rate can be approximated 
by a closed form expression presented in the following corollary. 

Corollary 1. (Closed Form Approximation of Individual Sum Rate) 

h^Ri™ = M±Q (-L)" (i-rp^Mr,). (22) 

Proof: The proof is given in Appendix |Bj ■ 
In order to demonstrate the rate approximation in Corollary [Tj we conduct a numerical study 
in Fig. [2] for different M and pu with respect to the number of users. The exact Rk in (1201) can 
be calculated by numerical integration. It is observed that (|22|) approximates the exact rate very 
well, which makes the rate approximation valuable due to its efficient computational form. 

B. Individual Scaling Laws 

The difficulty of dealing with rate scaling with selective feedback is two-fold. Firstly, due to 
selective feedback of the best beam, the number of SINR to maximize over at the scheduler side 
for each beam is a random quantity. This random effect is reflected in the random set U m in 
Section ITV-AI Secondly, the normalizing constants for establishing the type of convergence fl3TI . 
Il33l have to be obtained for a quantity d other than the number of users K in the full feedback 
case. In [|22l. the first issue was tackled by the Delta method. In this paper, we solve the first 
issue by referring to the extremes over random samples, and rigorously solve the second one by 
using the normalizing constants theorem. The proof is provided in Appendix El 

To examine the random effect on multiuser diversity, denote the sequence of random variables 
K m (K) as the number of SINR fed back for beam m with K users. It is easy to see that n m (K) 
are binomial distributed with probability of success -k. Thus by the strong law of large numbers, 
as K grows, the number of SINR fed back for each beam becomes jj. The following theorem 
is called upon to deal with this random effect. 
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Fig. 2. Comparison of the exact individual sum rate and the approximated one for a given user with different M and pt with 
respect to the number of users (M = 2,4, p k — dB, 10 dB, 20 dB). 



Theorem 3. {Extremes with Random Sample Size 071/ . Let, as K — v oo, — > i? in 

probability, where i9 is a positive random variable. Assume that there are sequences clk g 
M, 6^ > such that K \^ K converges weakly to a nondegenerate distribution function C7. Then, 

as K — > oo, 

/oo 
C7^(x)dP(tf < y). (23) 
-oo 

Therefore, if we denote Ak: K (K) as the extreme order statistics of the received SINR for each 
beam of a given user k, then from Theorem |3] its CDF can be efficiently approximated by (i*V fc ) 13 ■ 
Combining this with the normalizing constants theorem in Appendix |B] yields the following 
corollary. 
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Corollary 2. (Individual Scaling Laws Under Spatial Dimension Selective Feedback) 

p pApp 

lim tit = i i im n \ = 1. (24) 

K->oo M log 2 log 2 jj K^oo M log 2 log 2 K 

Proof: The proof is given in Appendix El ■ 
Remark: The scaling for the exact rate Rk and approximated rate -R^ pp differs in the factor 
■fij. The rate scaling for R^ pp does not have this factor because intuitively the exponent M in 
the approximated CDF (Fz k (x)) M counteracts the reduction in the number of SINR values for 
maximization, i.e., J|, due to selective feedback. We call this effect as the multiplicative effect. 
The detailed proof can be found in Appendix EJ To draw further insights, we can think of the 
exponent of Fz k (x) as the virtual users. In the full feedback case, the exponent equals K. In 
the selective feedback case with the approximated CDF, the exponent asymptotically equals K 
by the aforementioned multiplicative effect. The notion of virtual users and the multiplicative 
effect will be investigated further with both spatial and spectral dimension selective feedback in 
Section IV-Cl 

V. Selective Feedback in Both Spatial and Spectral Dimension 

In this section, random beamforming is embedded in a wideband OFDMA system. The system 
model is presented in Section IV-A[ the exact analysis and the asymptotic analysis are examined 
in Section IV-BI and Section IV-Cl respectively. 



A. System Model 

The system model described in Section [XT] is extended to an OFDMA system with N resource 
blocks. Each resource block is regarded as the basic scheduling and feedback unit. The random 
beamforming strategy generates M orthonormal beams 4> mn for each resource block. Denote 
s m n as the mth transmission symbol at resource block n, then the received signal y kjU for user 
k at resource block n can be expressed as 

M 

(25) 

771=1 



5 Note that even though the scaling laws are the same for the full feedback and the selective feedback case, this metric 
only measures the asymptotic performance when K is large. The exact rate performance is different due to the randomness of 
multiuser diversity and the scheduling outage event. 
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Fig. 3. Illustration of the spatial and spectral dimension selective feedback and the scheduling result in an OFDMA system 
(different colors denote different users K = 9, N — 5 resource blocks, M = 4 beams, the spectral dimension selective feedback 
L = 2). 



where H fc „ G C Mxl is the frequency domain channel transfer function of user k at resource block 
n with i.i.d. CA/"(0, 1) elements. To facilitate analysis, H fc n is assumed to be i.i.d. across resource 
blocks for a given user. This corresponds to the widely used block fading approximation in the 
frequency domain 11351 . [[36 J due to its simplicity and capability to provide a good approximation 



to actual physical channels. The transmit power for a resource block is assumed to be 1. From 

l H L*m,»| a 



(1251) . the SINR/c nm of user k at resource block n for beam m is SINR fc n m = — , — ; . 

and is denoted by Z kyU , m for notational simplicity. For a given user k, the Z\. n m 's are i.i.d. 
across resource blocks for a given beam m, and for a given resource block n, the Z k n >m 's are 
identically distributed and correlated across beams. The CDF of Z knm is given by F Zk (x) = 

-jr x \ 

1 — ) where the index n and m can be dropped due to the identically distributed 

property. 
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B. Individual Sum Rate 

With the extra degrees of freedom in the spectral dimension, additional selective feedback at 
each user side can be made possible by the following two-stage feedback selection. The first 
stage selection is in the spatial dimension, where each user selects the best beam with the largest 
SINR for each of the resource block. This process is similar to the narrowband feedback selection 
discussed in Section HY-Al Let Y^ n<m be the outcome of the first stage selection, thus from (fT5l ). 

( 2Mx \ 
1 — ^ X ^ + a{ x ) I u ( x )> where again the resource 

block index n and the beam index m can be dropped due to the identically distributed property 

across resource blocks and beams. The second stage selection occurs in the spectral dimension, 

where each user feeds back the SINR values of the best L resource blocks among the total 

N resource blocks. Let Wk, n ,m denote the outcome of the second stage selection of user k at 

resource block n for beam m. Thus this random variable represents the selected SINR at the 

user side, whose CDF is of interest for further analysis. It is easy to see that for the case of full 

feedback in the spectral dimension, i.e., L = N, F Wk = F Yk . For the best-1 feedback case, i.e., 

L = 1, Fw k = (Fy k ) N due to the independent property of Y k across resource blocks. For the 

general best-L feedback case, utilizing the results in [|29l , the CDF can be shown as 

L-l 

Fw k {x) = Y,UN,L,i)(F Yk (x)) N - e , (26) 
e=o 

where £i(N, L, £) = J2i=e ^JT (T) G) ( — 1)* ^ The two-stage feedback selection is demonstrated 
in Fig. [3] with nine users denoted by different colors, five resource blocks, and four beams. In 
the illustrated example, we use best-2 spectral dimension feedback, i.e., L = 2. 

After receiving feedback, the scheduler performs the CDF-based scheduling by first conducting 
the transformation on the received SINR, 

Wk, n ,rn — Fw k [Wk,n,rn) ■ (27) 

Denote k* n m as the random variable representing the selected user at resource block n for beam 
m, then 

k* nm = max W k , n ,m, (28) 

where W n>m denotes the set of users conveying feedback for beam m at resource block n. 
Following the derivation in Section |TV-A[ let X n>m be the selected SINR for beam m at resource 
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block n at the scheduler side. Then averaging over the randomness of \U n}7n \, the conditional 



CDF conditioned on k* n m can be written as 



A' 



ti=0 



K\ ( 1 



TiJ \M 



M 



K—t\ n 



E 



T 2 =0 



n\ ( l. \ ( l j N 

1 ~ N 



r 2 J \N 



(29) 

For further derivation, (Fw k (x)) T2 is manipulated into the following form by the power series 
expansion ||29l 11571: 



r 2 (L-l) 

(Fw k (x)r= UN,L,r 2 ,£)(F Yk (x)) N ^, 
e=o 



(30) 



where 



(£i(/V,L,0)r 2 , £ = 



&(N,L,t 2 ,£)= { 



(3D 



ttl(N,L,0) 

x&(N, L, i)&{N, L, t<2, £ — i), \<£<t 2 (L-1) 

y (d(N,L,L-i))^, e = r 2 (L-i). 

Following the same procedure as in Section IIV-A1 the individual sum rate for user k can be 
derived as 



1 - 
Rk = — y E 



n=l 
A 



71 = 1 

r 2 (L-l) 

X 

£=0 



M 

^log 2 (1 + X n>m \k* nim = k) 

m=l 

rJ I M 



A— 71 71 



UN,L,r 2 ,£) log 2 (l + x)d(Fy k (x)) 

v-n JO 



Nt 2 -£ 



(32) 



In order to obtain the closed form rate approximation for R k , the CDF approximation proposed 
in (|21"l) by the negative association property is utilized to approximate F Yk . The closed form 
result is presented in the following corollary. 

Corollary 3. (Closed Form Approximation of Individual Sum Rate) 



R k ~ R^ pp = M 



A / / 1 \ n / 1 \ A-Ti 71 

-) (i--) y 

M I \ Ml L^i 



71=1 

r 2 (L-l) 



£ &(N,L,r 2 ,£)J'k{M(NT 2 —£)) 



72=1 



Ti 

r 2 7 VAT 



/. \ ra / ( /. x 



(33) 



£=0 
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Fig. 4. Comparison of the exact individual sum rate and the approximated one for a given user with different spectral dimension 
selective feedback L with respect to the number of users (M = 4, N = 10, pk = 10 dB, L = 1, 2, 4, 10). 



To understand the impact of spectral dimension selective feedback, we conduct a numerical 
study assuming N = 10, M = 4. Fig. |4] plots the exact and approximated rate for different L 
under p k = 10 dB with respect to the number of users. It can be seen that when the number of 
users is small, there is a certain rate gap between selective feedback and full feedback. However, 
the gap becomes negligible when the number of users increases. In Fig. [51 the performance is 
observed for different p k for K = 20. From the two figures, we can see that the proposed rate 
approximation tracks the exact performance very well. 

C. Individual Scaling Laws 

We now examine the rate scaling with selective feedback in both spatial and spectral dimension. 
In Section IIV-BI with spatial dimension selective feedback, the CDF of interest is F Yk and the 
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Fig. 5. Comparison of the exact individual sum rate and the approximated one for a given user with different spectral dimension 
selective feedback L with respect to different p k (M = 4, N = 10, K = 20, L = 1, 2, 4, 10). 



number of SINR to maximize over at the scheduler side for each beam approaches j^. With 
additional spectral dimension feedback, the CDF of Fw k is of primary interest. To get a handle 
on the randomness of multiuser diversity for this case, an approach similar to that in Section 
IIV-BI can be utilized. Let the sequence of random variables K n ^ m (K) be the number of SINR 
values fed back for beam m at resource block n with K users. It is easy to see that n n>m (K) 
are binomial distributed with probability of success -r^. Therefore, by the strong law of large 
numbers, as K grows, the number of SINR values fed back for each beam at each resource block 
becomes -^4. Moreover, the convergence property of the sequence K n m (K) can be shown by 
invoking the central limit theorem: 
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Fig. 6. Comparison of the exact CDF Fw k and its tail equivalence for different spectral dimension selective feedback L 
(M = 4,N= 10, p k = 10 dB, L = 1,2, 4, 10). 



where d indicates convergence in distribution. By applying Theorem[3l the extreme order statistics 
of the received SINR for each beam at each resource block for a given user k can be efficiently 

KL 

approximated by (F Wk ) MN . 

Now the remaining problem is to examine the type of convergence of F Wk . Recall the 
formulation of F Wk as: F Wk (x) = Y^=q ^i{N,L,t)(F Yk {x)) N '^ . It is known that F Yk converges 
weakly to the type 3 Gumbel distribution. Due to the complicated form of •, •), it is tedious 
to directly check the conditions for proving the type of convergence. In order to investigate the 
tail behavior of F Wk which dominates the type of convergence ll33l . the following tail equivalence 
theorem is called upon. 

Theorem 4. {The Tail Equivalence Theorem /IMI/) U(-) and V(-) are distribution functions such 



IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS, TO APPEAR 



21 



that 

lim = L (35) 

x-+oo 1 — V (X) 

If there exist normalizing constants a K , bx > such that U K (a K + b^x) — > G(x), where G(x) 
is non-degenerate, then V K (aK + bxx) — > G(x). 

From Theorem 0] one can infer that if two distribution functions are tail equivalent, then they 
belong to the domain of attraction of the same type. Employing Theorem HI a tail equivalent 
formulation can be obtained for F Wk expressed in the following corollary. 

Corollary 4. {Tail Equivalent CDF) F Wk (x) is tail equivalent to (F Yk (x)) N -^e=o fiW^. 

Proof: The proof is given in Appendix O ■ 
Fig. |6] compares the exact CDF and the corresponding tail equivalence for different selective 
feedback L under M = 4, N = 10, and pk = 10 dB. The tail equivalent CDF is observed to track 
the exact one even when x is small, which supports and lends confidence in the power of the tail 
equivalence theorem. Therefore, the tail equivalence is used to study the type of convergence, 
which is expressed in the following lemma. 

Lemma 2. (Type of Convergence of Selective Feedback) Given the statistical property of F Yk in 
(fl"5l) . Fy/ k belongs to the domain of attraction of type 3 Gumbel distribution. 

Proof: The proof is given in Appendix O ■ 
Having obtained the type of convergence for Fw k , the rate scaling result can be derived by 
referring to the normalizing constants theorem in Appendix [Bj The individual rate scaling is 
provided below. 

Theorem 5. (Individual Scaling Laws Under Spatial and Spectral Dimension Selective Feedback) 

Hm Rk = 1 Hm b£™ = 1 

K ^ M log 2 log 2 W-tttf^ K ' ' *™ M lQg2 lQg2 (N-^ ( U,(N,L, m L K 

(36) 

Proof: The proof is given in Appendix O ■ 



Remark: For the exact rate Rk, the ultimate equivalent CDF of interest is F Y 



K 



( N— y^ L ~ 1 £i (N L £)£)L 

thus the exponent - — a/at ' K due to multiplicative effect can be seen as the virtual 
users for scheduling competition. This exponent is for the general best-L spectral dimension 
feedback. For the full feedback L = N, since £i(N, N, €) equals 1 for £ = N — l and otherwise, 
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K 



the CDF becomes Fy . For the best-1 feedback L = 1, since ^(N, 1, £) is 1 for i = and 



otherwise, the CDF becomes Fy* . Intuitively, the best-1 feedback is asymptotically optimal 
due to the same number of virtual users. In other words, even though additional maximization 
reduces the average number of variables for feedback, it counteracts this reduction by increasing 
the exponent of the CDF The number of virtual users is the limiting factor that dominates rate 
scaling. For the approximated rate R^ pp , since the approximated CDF compensates for the spatial 
dimension selection by increasing the exponent, the rate scaling differs by a factor of M. 



In this paper, an analytical approach is used to investigate the problem of random beamforming 
with heterogeneous users and selective feedback. The heterogenous user scenario corresponds to 
the practical scenario of potentially different large scale channel effects for different users. We 
leverage the CDF-based scheduling policy to decouple each user's rate and thus theoretically 
examine the individual user rate. We develop the notion of individual sum rate to analyze the rate 
scaling for each individual user. We focus our analysis in this work on theoretically understanding 
the effect of selective feedback in both spatial and spectral dimensions. On the exact analysis 
part, extensive numerical results show that our approximate expression for the rate under selective 
feedback is effective and provides an efficient expression for computing the exact rate. On the 
asymptotic analysis part, we develop the notion of virtual users and the multiplicative effect to 
explain the impact of selective feedback on rate scaling. We further discover that the limiting 
factor for the rate scaling is the exponent for the ultimate CDF of the selected SINR at the 
scheduler side. The extension of this work from single antenna users to multiple antenna users 
and more generally the investigation of a multicell MIMO setup with advanced receiver design 
will be the subject of our future work. 



K 



VI. Conclusion 



Appendix A 



Proof of TheoremUl' With the help of Lemma [H J7fe(e) can be computed as 




) 



(37) 
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where (a) follows from applying Lemma [IJ (b) follows from integration by parts. The closed form 
result for Z(a, /3) in Theorem Q] can be computed in a recursive manner [|3~7l and is presented 
as follows 



Z(a,f3)={ (/3 " 1)! ti^- ' (38) 



e a E 1 (a), (3 = 1 

where E 1 (x) = ^-r-dt is the exponential integral function of the first order []39 

Appendix B 

Proof of Corollary [7} 



Mr, 



where (a) follows from the CDF approximation in (l2TT ); (b) follows from the definition and 
computation of J^ie). 

Proof of Corollary\2\- It is shown in [22] that Fy fc belongs to the domain of attraction of type 
3 Gumbel distribution ll33l . Thus if the number of SINR to maximize over for each beam is 
fixed and equals the number of users K, then the following equation holds: lim (F Yk (ak-.K + 

K—>oo 

i>k-.Kx)) K = &(x), where ty(x) = e~ e * is the type 3 Gumbel distribution, a^-.K and b k:K represent 
the normalizing constants for user k. From Theorem |3] the number of SINR to maximize over for 
each beam approaches ||. Let c^-.k an d dk,K denote the normalizing constants for user k under the 



selective feedback case. Then the following equation holds: lim (FyJck-.K + dk-.Kx)) M = ^f(x). 

K—>oo 

In order to obtain c^-.k and dk,K, the following theorem is called upon. 

Theorem 6. (The Normalizing Constants Theorem 071/) Let Fk{v) be a sequence of distribution 
functions. Let a K , b K > 0, ck, and d K > be sequences of real numbers such that, as K — >■ oo, 

\mxF K {a K + b K x) = U(x), lim F K (c K + d K x) = V(x) (40) 

for all continuity points x of the limits, where U{x) and V(x) are nondegenerate distribution 
functions. Then, as K — > oo, the limits: lim^ = B ^ 0, lim CK ^ K = A are finite, and 
V(x) = U(A + Bx). 
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The spatial dimension selective feedback case possesses the following situation in Theorem 

[6l F K (x) = (F Yk (x)) K , a K = a k:K , b K = b k:K , c K = c k]K , d K = d k:K , U(x) = $(x), and 
V(x) = ($f(x)) M . The sequence of a k:K has been derived in [12] as: a k , K = p k log 2 K — p k (M — 
1) log 2 log 2 K + o(l). A suitable choice of b k , K for type 3 is g k (b k .x), where g k (x) is the growth 
function for user k defined by g k (x) = 1 f z f ^ . Thus a suitable sequence is b k: K = pk for all K. 

JZ k \ x ) 

Solving (ty{x)) M = ty(A + Bx) yields A = — log M, B — 1. Therefore, by referring to Theorem 
(6[ the normalizing constants can be derived to be: c k:K = p k log 2 J|— p k (M— 1) log 2 log 2 K+o(l), 
and d k:K = p k for all K. Then by employing the Corollary A.l. in [12], the individual rate for 
user k, namely R k scales as Mlog 2 log 2 

Regarding the approximated rate R k pp , since the approximated CDF by negative association is 



(F Zk (x)) and the number of SINR to maximize over approaches 4|, we have lim (F Zk (c k: K + 
d k: Kx)) MjI = lim (F Zk (c k ;K + d k: Kx)) K = *f?(x). Thus the normalizing constants c k: K = a-k-.K, 

K—>oc 

and d k: K = b k ._K, which enables the approximated rate R^ pp to scale as M log 2 log 2 K. 

Appendix C 

Proof of Corollary^ Given Theorem @] the following equality holds: 

lim l-EUm,L,l)(F Yk (x)r-t 

x->oo 1 - (Fy (x)) N -^=o £i(JW)< 



(41) 



(A r - Et"o 1 ei(^,L,£)£)(Fy fc (x))^-fe 1 6(^K-i /yfc(x ) ' 

where (a) follows from the L'Hospital's rule; (b) follows from the fact that J2e=o £i(-^> ^) = 1- 

Proof of Lemma \2\- F Yk with statistics in (T15T ) belongs to the domain of attraction of type 
3. It can be shown that for any distribution function F(x) which converges weakly to the 
limiting distribution, then its exponent form F £ (x) has the same type of convergence [1311 . 
(F Yk (x)) N ~^e'=o ^( N ' L ' e ) e belongs to the domain of attraction of type 3. Then by Theorem |4] 
F Wk belongs to the domain of attraction of type 3. 

Proof of Theorem\5\ A procedure similar to that used in proving Corollary 0can be used here. 
Since the number of SINR to maximize over for each beam at each resource block approaches 



, and F Wk belongs to the domain of attraction of type 3, the following equation holds: 



KL 

MN 

lim (F Wk (c k:K + d k ; K x))^f = ty(x). By referring to the tail equivalence theorem, the equivalent 

K— >oo 
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equation is: lim (F Yk (c k:K + d k . K x)) mn = v]/(x). Applying Theorem [6] yields 

K—¥oo ' 

the normalizing constants: c k , K = p k log 2 KL ^ N ~Ei=o^ N ' L ' t > t > _ p fc (M — 1) log 2 log 2 K + o(l), 
and d k ., K = p k for all K. Therefore, R k scales as M log 2 log 2 KL( - N ^^^ N ' L,t>i \ 

For the approximated rate R^ pp using the approximated CDF (F Zk {x)) M for F Yk , the following 
equation holds: lim (F Zk (c k: K + d k -Kx)) ^ = Using the same line of 

K— >oo 

arguments, i?£ pp scales as Mlog 2 log 2 W^E^M^M , 
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